Intrinsic Tunneling Spectra of Bi 2 Sr 2 CaCu 2 08+5 near Optimal Doping 
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We report tunneling spectra of near optimally doped Bi2Sr2CaCu208+<5 intrinsic Josephson junc- 
tions with area of 0.09 /im 2 , which avoid some fundamental difficulties in the previous tunneling 
experiments and allow a stable temperature-dependent measurement. A d-wave Eliashberg analysis 
shows that the spectrum at 4.2 K can be well fitted by considering electron couplings to a bosonic 
magnetic resonance mode and a broad high-energy continuum. Above T c , the spectra show a clear 
pseudogap that persists up to 230 K, and a crossover can be seen indicating two different pseudogap 
phases existing above T c . The intrinsic electron tunneling nature is discussed in the analysis. 

PACS numbers: 74.50.+r, 74.25.Jb, 74.72. Hs 
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The paring mechanism and pseudogap phenomena are 
known as the two key issues in the studies of high-T c 
superconductors. In addition to the views based on res- 
onating valence bond (RVB) state [lj], one possible mech- 
anism that is of much interest is boson-mediated pair- 
ing of the electrons. Recent angle-resolved photoemis- 
sion spectroscopy experiments have focused on a "kink" 
in single-particle dispersion and both phonon 0, Q and 
magnetic resonance [J [B| are considered as the bosonic 
candidate. In the tunneling and optical Q experi- 
ments, a broad high-energy continuum, in addition to 
the resonance mode, has also been considered. 

Tunneling has traditionally been a useful tool in re- 
vealing the material's superconducting properties. For 
the high-T c superconductors, mostly Bi2Sr2CaCu20s+5 
(Bi-2212), the scanning tunneling microscope (STM) [8|, 
break junctions (BJs) la lSjlK and intrinsic Josephson junc- 
tions (IJJs) [13, m Hill 
many important informations 
ski et al. have correlated the tunneling dip structure 
with the magnetic resonance mode @. In Ref. [H, the 
pseudogap opening temperature T*, rather than the su- 
percondcting transition temperature T c , is suggested to 
be the mean-field critical temperature, thus supporting 
precursor pairing as the origin of the pseudogap phase. 

IJJs, being intrinsic, have the advantage of avoiding 
the problems like the surface chemical deterioration and 
unstable junction structure, and can offer a convenient 
temperature-dependent measurement. In the previous 
studies, however, the micro-meter-sized IJJs are found 
to suffer from significant self-heating, which severely dis- 
torts the tunneling spectra H, 14 1 . 

Recently, we have demonstrated that self-heating in 
IJJs can be reduced considerably when junction sizes de- 
crease down to the sub-micron level [14j. In this Letter, 
we report and present a detailed analysis on the tunnel- 
ing spectra of an IJJ mesa with further-reduced area S 
= 0.3x0.3 /j,m 2 on a near optimally doped Bi-2212 crys- 
tal which are shown in Fig. [T] Comparison with 
the previous reports (including those from STM and BJ 
experiments) is presented in Table [I] in two important 



spectroscopic aspects: the relation between the super- 
conducting gap (SG) and pseudogap (PG), and whether 
there is a clear dip at low temperatures. It can be seen 
that self-heating can lead to considerably different results 
in the IJJ experiments. 

The spectra in Fig. [1] show a well-defined, almost 
temperature-independent normal-state resistance Rn, 
which allow a straightforward tunneling analysis. Below 
we analyze the data, along with a discussion on the tun- 
neling nature in IJJs, of which a consensus is still lacking. 
We will show that the spectrum at 4.2 K can be explained 
within the d-wave Eliashberg description by considering 




FIG. 1: (color online). Temperature dependence of the I-V, 
dl/dV, and d 2 I/dV 2 curves of a Bi-2212 mesa (T c = 89 K) 
containing 10 IJJs. V corresponds to the voltage per IJJ. 
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TABLE I: Some differences and similarities of the tunneling 
spectra for (near) optimally doped Bi-2212. Results for STM, 
BJ, and IJJ are from Refs. @, and [Til, respectively. 





STM 


BJ 


IJJ 


This work 


SG vs. PG a 


SG^PG 


no PG 


SG/PG 


SG^PG 


Clear dip? 


yes 


yes 


no 


yes 



"SG^PG and SG/PG denote a smooth evolution from SG to PG 
at T c and that they coexist and evolve independently, respectively. 



a magnetic resonance mode plus a broad continuum [6 1 , 
with the latter having an upper cutoff at 300 meV [17|. 
Above T c , there is a clear indication of the existence of 
two pseudogap phases, which has not been reported so 
far in the tunneling experiments for the (near) optimally 
doped Bi-2212 samples. 

The d-wave Eliashberg formalism. — We use a boson- 
exchange spectral density in a separable form (l8j : 



*k,k 



F(Q) =a 2 [c s + c D cos(29)cos(29') }T(fl) , 



where k or 9 and k' or 9' are associated with the initial 
and final states in pairing interaction, c$ and crj denote 
the s- and d-wave com pon ents, and J-(Cl) is the cut-off 
Lorentzians defined by [19( 
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\fi-n \ < fii , 



\fi-n \ > fii , 



in which C normalizes !F(Q) to unity. The Eliashberg 
equations in the d-wave form at T — can be written as: 



oc 2tt 




oc It: 



A(uj')cos 2 29'K + (uj,u>') 
^Juj' 2 - A 2 (uj')cos 2 29' 



V^lu' 2 - A 2 (uj>)cos 2 29> 



in which A and Z are the gap and renormalization func- 
tions, and K±(u),u>') is given by 

oo 

K±(uj,uj') = [ dfl a 2 T(Q) ( —— 
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The above equations can be solved if the parameters cs, 
crj, and a 2 J-(VL) are known. Once A(w) and Z(uS) are 
obtained, the diagonal self-energy of the system can be 
found from ReE = u>[l— ReZ] and Im£ = — ojlmZ, and 
the density of states (DOS) is given by 

N(u3,6) = Re[u/^u 2 - A 2 (uj)cos 2 (29)} . 



Discussion on the tunneling nature in IJJs. — A usual 
way to fit the tunneling spectra is to build an angle- 

averaged DOS: N d (uj) = (l/2vr) / N{uj,9)d9, from which 

o 

the I(V) curve can be calculated straightforwardly con- 
sidering an incoherent tunneling process @, 0, 21 1 . The 



basic features from such calculation at 4.2 K using a set 
of parameters described below can be seen in Fig. [2a) 
and (b) (dotted lines). While the lines follow the exper- 
imental data closely above the gap voltage, they devi- 
ate considerably below it. From our calculation, we find 
that a better fit below the gap voltage can be obtained 
by considering coherent electron tunneling with a direc- 
tional matrix element of t ~ cos 2 (29), but this leads to 
a worse agreement near and above the gap voltage [2fj| . 

Coherent electron tunne ling has been discussed in a 
number of publications [l^, 22, 23 1 and used to interpret 
the c-axis transport data in IJJs 24, (2||. The present 
results may imply that there is a crossover from coher- 
ent to incoherent tunneling as the bias voltage increases. 
Difficulty arises, however, that when a ( ~ sin 2 (49) is 
used, which is theoretically more strict for the Bi-2212 
material in which the Cu atoms in Cu02 planes across 
the BiO and SrO barrier layers do not lie collinearly [23| , 
the computed results give a worse fit for voltages both 
below and above the gap edge. To avoid the concep- 
tual difficulty, we assume that intrinsic electron tunnel- 
ing in IJJs is basically incoherent for voltages away from 
zero. At low bias, tunneling can be described by an effec- 

tive DOS: N eff (w) = (1/tt) J N(u,9)cos 2 (29)d9, which 

o 

produces similar I-V and dl/dV characteristics (dashed 
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FIG. 2: (color online). Fit of the experimental (a) I-V and 
(b) dl /dV curves at 4.2 K using the d-wave Eliashberg theory. 
The calculated DOS and self-energy E are shown in (c) and 
(d). Note the horizontal shifts of 150 and 300 meV for N R ff 
and N(uj,0) in (c) for clarity, a 2 T used in the calculation is 
also shown in (d). See text for detailed fitting parameters. 
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lines in Fig. Ufa) and (b)) as coherent tunneling plus t ~ 
cos 2 (29) is considered. Physically, the cos 2 (26) factor in 
N e ff can result from the wave-function overlap (having 
d x 2_ y 2 symmetry) between the 0-2p and Cu-4s orbitals 
while the latter plays an important role in assisting a 
hole hopping to and out of the 0-2p orbital [12] ■ At high 
bias, tunneling is described by Nd(uj). 

Switch from N e f / to Nd implies a lesser role the Cu-4s 
orbital plays in assisting hopping as bias increases and 
more high-energy quasiparticles with reduced lifetimes 
are involved in tunneling. Phenomenologically, the tun- 
neling current can be written as I = PId + (1 — P)I e ff, 
where Id and I e f f represent those obtained using Nd and 
N e ff, and P is a voltage-dependent probability distri- 
bution. In our fitting, an exponential function P(V) — 
(e v / Vl -l)/(e V2 / Vl -l) with < V < V 2 is used, in which 
V\ characterizes the voltage scale for I e f f to Id crossover 
and Vi sets a point above which I = Id- The solid lines 
in Fig. [21(a) and (b) are obtained from the dotted and 
dashed results with V\ and V2 equal to 20 and 75 mV, 
respectively. The fit appears satisfactory. 

Tunneling spectrum at 4-2 K. — In the Eliashberg anal- 
ysis of the 4.2-K data in Fig. [21 we follow Zasadzinski et 
al. and use two Lorentzians for a 2 J- to model the mag- 
netic resonance peak (around 35~40 meV) and the high- 
energy continuum The fitting parameters are: cs — 
0.18, c D = l-a\ = 20, fi i = Oil = 35 mcV, 771 = 5 meV; 
a\ = 46.5, fi 02 = fii2 = 150 meV, 772 = 2000 mcV; R N 
= 1.92 kSl. The parameters are adjusted so as to have a 
fair fit of the voltage positions of the dl/dV peak and dip 
and also the depth of the dip, which are therefore judged 
primarily by Id- The resulting antinodal DOS N(u>,0), 
together with Nd and N e ff, is shown in Fig. Etc), and 
o?T is depicted in (d). 

o?T is found to change the shape of the dl/dV curve in 
the following way. While the resonance peak creates the 
dip, the high-energy continuum tends to reduce its depth. 
In addition, as the continuum extends to higher energy, 
the width of the dip decreases. The latter change is more 
obvious when the upper cutoff of the continuum reduces 
from 150 to 80 meV, where the resulting dip width in 
both the dl/dV and DOS curves becomes too wide to 
fit the experimental data. This means that a continuum 
extending above 150 meV must be considered. 

We have chosen an upper cutoff of 300 meV for the 
continuum following a recent analysis of the optical data 
17| . As can be seen in Fig.[2td), this results in a consider- 
able high-energy contribution to — 2Re£, which may need 
further experimental testing. Using a smaller cutoff of ~ 
160 meV [6j will reduce such contribution and slightly in- 
crease the dip width. The overall size of — 2Re£, on the 
other hand, is mostly affected by cs- Reducing cs will 
reduce — 2ReS and make the peak-dip separation smaller. 

Compared with the BJ results [f|, the present spectra 
show a shallower dip, which leads to a larger contribution 
of the continuum in a 2 T relative to the resonance mode. 
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FIG. 3: (color online). Experimental dl/dV curves near and 
above T c (symbols). Lines are fits using N e ff, in which 
N(u,6) = Re[(u - iF)/y/(uj - iT) 2 - A 2 cos 2 (26)}, with A 
and r values plotted in Fig. [4] Curves are shifted downwards 
successively by 0.3 fcSl -1 for clarity except for that at 75 K. 



Presently, it seems unclear as to which part in o^T should 
be more essential for the pairing @, 0] . 

Spectra near and above T c . - As temperature in- 
creases, the dip gradually disappears at ~ 75 K (see 
Fig. QIb) and (c)). In this case, it is possible to fit the 
data using N e ff in which N(ui,&) has a Dynes' form 
Re[(u - iT)/y/{u - iT) 2 - A 2 cos 2 (20)} where A and T 
are the energy-independent parameters [26| . The fitting 
can be satisfactory using N e f f in the entire voltage range 
since difference of using N e f f and Nd above the gap volt- 
age is smaller at higher temperatures and in the absence 
of the dip. Comparison between experiment and theory 
can be seen in Fig. [31 and the A and T parameters are 
plotted in Fig. [U It is interesting to note the linear tem- 
perature dependence of the scattering rate T, a signature 
predicted in the marginal Fermi- liquid theory [27] . 

It should be pointed out that using A to describe 
the data above T c implies that we are taking the pre- 
pairing view and A is understood as the thermodynamic- 
averaged value j28|. The pre-pairing idea has been sug- 
gested previously in the STM experiments to account 
for the pseudogap phase [§]. From our experiments, the 
dl/dV spectra show a gapped structure for temperatures 
as high as 230 K. To look at this in more detail, we plot 
in Fig. 0] the temperature-dependence of half the peak 
energy in the dl/dV curves (Ep/2: solid squares). It 
can be seen that Ep/2 decreases as the temperature in- 
creases up to T c . It then increases, reaches a maximum 
at ~ 140 K, and starts decreasing again after that. The 
open squares in the figure are the results taken from the 
dl/dV curves that are normalized to the 250-K curve. In 
this temperature decreases from 230 K, Ep/2 is 
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FIG. 4: (color online). Peak energy Ep divided by two from 
experimental dl /dV curves (solid squares). Open squares are 
the results from the dl/dV curves normalized to the 250-K 
curve. Solid and open triangles are the A and V values used 
for fitting in Fig. |1 The dashed line is A(T) from the BCS 
d-wave description with 2A(0)/fcT c = 4.3. The dotted line is 
the corresponding result with temperature normalized to T c . 

almost flat until ~ 170 K, below which it decreases at a 
relatively fast speed. 

These results point to a change within the pseudogap 
phase, which can be characterized by two temperatures 
T* = 230 K and T 2 * = 140-170 K. Considering our near 
optimally doped sample, this is consistent in the phase di- 
agram with the recent Nernst experiment [29j . From the 
RVB scenario [l[ , the results can be viewed as the opening 
of a spin gap below T*. Below , holons start to Bose- 
condense and the Nernst state results. A global supercon- 
ducting state is finally developed below T c through the 
Berezinskii-Kosterlitz-Thouless transition. In the view 
of boson-mediated pairing, the change around TJ should 
correspond to the point below which the superconduct- 
ing order starts to form out of (with T comparable to A 
at initial stage), and compete with the spin-gap state. 

Summary. - Tunneling spectra of near optimally 
doped Bi-2212 IJJs in the temperature range of 4.2 to 250 
K are presented. We have analyzed intrinsic tunneling 
process and used (i-wave Eliashberg theory and Dynes' 
DOS to fit the tunneling data. The fit between theory 
and experiment is satisfactory. We have shown that the 
electronic coupling to a magnetic resonance mode plus a 
broad continuum can well explain the spectrum at 4.2 K, 
thus providing an evidence for the electronic-originated 
boson- mediated pairing. However, a doubt remains as to 
which part in the boson spectrum is more essential for 
superconductivity. We have also found that there exist 
two pseudogap phases above T c . 
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